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Lagrangian  Boundary  Equation  and  Unsteady  Separation  of 
Elliptical  Cylinder  at  Angle  of  Attack 
Chen  Yunming 

(Institute  of  Mechanics,  Academia  Sinica) 


Abstract 


A  4a*5V  V  ' 


J 


this  work,  the  Lagrangian  boundary  equation 


is  introduced  to 


Steady  the  separation  of  an  abruptly  moved  elliptical  cylinder  at  the 
angle  of  attack.  Not  only  the  occurrence  of  separation  singular 
points  but  also  their  forward  movement  with  time  were  calculated. 
The  computation  also  revealed  that  the  separation  pattern  would 
change  from  trailing  edge  separation  to  leading  edge  separation. 
In  regard  to  the  methodology  of  computation,  a  transformation  was 
introduced  to  map  the  non-uniform  mesh  on  the  physical  plane) to 
the  uniform  mesh  on  an  auxiliary  plane.  Thus,  it  is  only 
required  to  tighten  the  mesh  (locally  on  the  physical  plane  to 
improve  the  accuracy  without  increasing  too  much  computing  time. 
In  addition,  the  advantages  of  the  equal  space  finite  difference, 
such  as  simplicity  and  high  accuracy,  will  remain.  Furthermore, 

'a  complicated  variable  boundary  problem  (separation  point  moving 
forward)  was  treated  as  a  fixed  boundary  problem  which 

j 

significantly  simplifies  the  process.  •  -*►  f  *  '  ■  ' 


I.  Boundary  Layer  Equation  in  Lagrangian  Coordinates 
In  fluid  dynamics,  the  Euler  descriptive  method  is  most 
commonly  used  because  the  majority  of  meaningful  flow  is  steady. 
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Furthermore,  the  stress  tension  of  a  Newtonian  fluid  in  Euler 
coordinates  is  simple.  However,  the  Lagrangian  method  is  also  very 
useful  for  problems  concerning  separation  with  unsteady  boundary 
layer,  non-Newtonian  fluids  whose  viscous  tension  is  related  to 
the  microscopic  history,  or  when  specific  particles  in  the  fluid  need 
to  be  tracked  individually  (such  as  turbulent  diffusion). 

In  order  not  to  complicate  the  problem,  the  discussion  is 
limited  to  the  unsteady  two-dimensional  boundary  layer  of  an 
incompressible  fluid.  The  following  basic  relations  can  be 
obtained  from  *?VVX-1 ,  Vx+Vx=0' 

V.  “  —yilJ  V,  —  rtIJ  ( 

where  Jsx^y^-x^y^ .  In  this  case,  the  continuity  equation  is 

+  ‘'i’N  “  0  (A) 

From  equation  (1)  and  u=xt v=yt ,  we  get 

4-  ('■>•.  -  w  -  If  -  •  <B) 


Thus,  the  continuity  equation  becomes 

J  =  x^yn  -  x^y^  =  1  (2) 

Because  the  Jacobi  determinant  represents  the  area  ratio  of 
mapping,  the  above  formula  is  the  Lagrangian  expression  of  the 
conservation  of  fluid  volume. 

Note  that  =  x^  -x^-—-  =  jj—  then  the  boundary  layer 

equation  becomes 


Manuscript  received  on  April  5,  1983 

1)  This  paper  is  dedicated  to  the  15th  anniversary  of  the  death 
of  Professor  Guo  Yonghui. 


We  get  the  following  by  expanding  the  parenthesis 


- — />.  ■+■  vlx!(U^r  —  +  x,utt  1 

p 


O') 


Note  that  equation  (3)  is  not  coupled  to  equation  (2).  It  can  be 
used  independently  to  solve  x(£,n»t).  In  this  case  equation  (2) 
can  be  used  to  find  y(C,n,t).  The  characteristic  functions  for 
equation  (2)  are 


dj 

dl 


di 


dy  - 

**’  TT-1’ 


(A) 


Hence,  we  know  that  dx/dx=0.  Thus,  we  can  integrate  along  the 
x^constant  curve  to  get 

-  [*’*  -  r*”’  (4) 

■  y  Jfi®-,,  Je*  x, 

2  2 

Let  us  assume  that  x=n  9(5)  and  5=(5+U0t)/n  ;  it  is  then 
possible  to  obtain  the  similarity  solution  to  equation  (3’)-the 
steady  motion  of  a  semi-infinite  plate  (Blasius  solution).  From 
equations  (2)  and  (3)  it  is  also  possible  to  get  the  Falker-Skan 
solution  and  the  merging  point  flow.  Nevertheless,  the  most 
noticeable  advantage  of  the  Lagrangian  coordinate  is  in  the 
treatment  of  unsteady  boundary  layer  separation.  With  the 
exception  of  the  classical  case  of  two-dimensional  steady 
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boundary,  the  criteria  for  boundary  layer  separation  are  still  not  yet 

clear.  For  instance,  there  are  many  types  of  separation  of  a  three- 

r  1  ] 

dimensional  steady  boundary  layer  .  There  were  controversies 

on  the  separation  of  an  unsteady  two-dimensional  boundary  layer. 

In  the  late  fifties,  based  on  special  cases  (separation  of  steady 

boundary  layer  in  solid  wall  motion),  Sears  et  al  introduced  the 

[21 

M-R-S  condition:  u=0  and  8u/0y=O  at  separated  points  .  In  the 
seventies,  Sears  and  Tellionis  used  the  divergence  of  numerical 
methods  to  determine  the  occurrence  of  singular  points- the 
physical  appearance  of  separation1  *  .  For  instance,  at  a 

certain  time  after  the  cylinder  begins  to  move  (t=2u  t/R=1.3)  the 
calculation  of  a  certain  time  in  the  flow  field  is  divergent,  a 
singular  point  appears.  However,  the  computation  performed  by 
Cebeci  et  al  could  last  until  t=2.8^^.  They  believed  that  the 
divergence  of  the  computation  performed  by  Tellionis  et  al  is  a 
problem  of  the  specific  scheme  and  does  not  necessarily  have  any 
physical  significance.  The  problem  became  clearer  after  Shen 
Shenfu  et  al  introduced  Lagrangian  coordinates  to  treat  the 
problem*^ . 

In  Lagrangian  coordinates,  a  steady  flow  is  also  unsteady. 

Thus,  it  is  possible  to  use  one  description  for  both  steady  and 
unsteady  boundary  layer  separations.  From  equation  (4)  one  can 
see  that  when  a  certain  point  in  the  flow  field  has  the  following 
properties 

x=  0 ,  x  =  0  (5) 

^  n 

then  the  integral  along  the  x=xg  curve  through  this  point  is 
going  to  approach  infinity  (on  the  1 l\f  Re  scale).  The  above 
formula  shows  that  the  fluid  masses  near  (5  ,n  )  reach  the  same 
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position  x*xg  in  the  x-direction  at  the  same  time.  Because  the 
volume  is  invariant,  thus  it  extends  to  infinity  in  the  y- 
direction.  It  seems  to  have  an  unsur mountable  obstacle  in  the  x- 
direction.  The  fluid  masses  cannot  continue  forward.  They  are 
separated  from  the  body  and  move  toward  infinity  in  the  y- 
direction  (on  the  scale  of  1/v/Ei).  This  is  the  definition  of 
separation  given  by  Prandtl  earlier. 

In  the  following  we  will  prove  that  equation  (5)  is  valid 
for  separation  of  a  steady  flow  in  a  solid  wall  movement.  It  is 
obvious  we  have 


Jtu(x,  4>) 


t(x , 


<l>,  S) 


(6) 


/,!{  H"  /,</<;  +  *{  —  0 

t, *,  +  ™  o 


(7) 


If  u(x  ,<p )=0  in  the  flow  field  (not  in  the  solid  wall  x4=0, 
s  s 

from  equation  (6)  we  know  that  t  -•  ®.  Based  on  equation  (7), 

equation  (5)  must  be  satisfied.  We  also  have  u  =x_u  -x  u  =0, 

y  5  n  n  5 

i.e.,  the  M-R-S  condition. 

In  Lagrangian  coordinates,  the  separation  condition  (5) 
shows  that  a  separation  point  is  an  extremum  of  x(g,n),  instead 
of  a  singular  point  (a  singular  point  appears  only  when  we  try  to 
determine  y  and  y  is  not  needed  in  calculating  x).  Therefore, 
there  is  no  difficulty  and  ambiguity  in  the  numerical 
computation.  It  is  clearly  indicated  in  reference  [6]  that 
cylinder  separation  begins  to  emerge  at  t=3.0  and  0g^111°. 
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II.  Unsteady  Separation  of  Elliptical  Cylinder  with  Angle 
of  Attack 

In  the  following  we  will  extend  the  above  one  step  further 
to  apply  the  Lagrangian  coordinate  to  calculate  the  unsteady 
separation  of  an  abruptly  moved  elliptical  cylinder  with  angle  of 
attack.  There  are  two  objectives:  first  we  not  only  want  to 
calculate  the  occurrence  of  an  unsteady  separation  but  also  to 
track  it  with  time.  Next,  we  want  to  study  the  variation  of  the 
position  of  separation  with  angle  of  attack.  That  means  when  the 
angle  of  attack  is  increased,  it  changes  from  trailing  edge 
separation  to  leading  edge  separation  (which  is  the  cause  of 
stalling  of  aircrafts).  As  the  first  step,  we  temporarily 
neglected  to  reverse  influence  of  separation  on  the  pressure 
distribution,  which  is  allowed  in  the  initial  stage  of 
separation.  This  reverse  influence  is  taken  into  account  in  the 
next  step.  Because  the  calculation  of  pressure  distribution  is 
an  independent  subroutine,  taking  the  reverse  influence  into 
account  has  no  effect  on  the  numerical  solution  of  the  boundary 
layer  itself.  It  is  essentially  a  problem  to  accurately  simulate 
the  separation  vortex  in  the  external  flow,  which  is  another 
special  topic  to  be  specially  discussed  elsewhere. 

In  the  following,  physical  quantities  with  '  represent 
dimensional  parameters  and  those  without  '  are  dimensionless 
quantities.  First,  we  introduced  an  elliptical  coordinate  to  use  <p 
to  replace  the  arc  length  x'  to  describe  the  points  on  the 
surface  of  the  ellipse,  i.e.,  to  change  variables  from  (5,n)  to 
(<*„  ,n).  The  unknown  function  is  still  x  (which  is  a  function 


d— St  u  j<i 
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of  9).  In  a  Cartesan  coordinate  (x,y),  an  ellipse  is  x=a  cos<p  and 
y*b  sin»  (where  a  and  b  are  the  long  and  short  axes).  The  arc 
length  of  the  ellipse  is 


■*  j - 

t'  —  V  a!sinJ</i  ■+■  i’cos’c/i  dtp, 

_  0 


It  is  rendered  non-dimensional  as 


•—  — — —  “  f  \/ E„  +  E2  sin" 
c j  -f-  b  J  o 


where  Eo=(2b/a+b),  and  E2  =4  a-b/a+b.  The  inverse  is  noted  as  9 
(x).  In  this  case  the  variables  are  £=x(<p0  )  and  90=p(5).  The 
non-dimensional  potential  velocity  is 


U (  4 1)  — •  u  —  cos <t  sin  </«  sin  cos 

■%/ E,,  ■+■  E,  slu'd/ 


where  a  is  the  angle  of  attack.  The  origin  is  placed  at  the  apex 
of  the  a  semi-axes.  Therefore,  the  leading  edge  stagnation  point 
is  9=-o,  x# =x(-d)=-x(a) .  The  time  is  rendered  dimensionless  as 
t=4u  t/a+b.  The  ordinate  is  also  made  non-dimensional  as 

09 

y=2y'/a+b  \Zua(a+b)  /  v  .  The  non-dimensional  velocity  is 
u=u ' /2u  . 

OP 

Next,  the  following  transformation  is  introduced  to  map  the 
semi-infinite  plane  (-d^<p0  <  *— o,0^n<+® )  to  the  square  (O^a^l,  0^8 
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(10) 


where  C,  hj,  and  h2  are  constants.  hx  =4(  1+d) /  (srA,(2  +  <0)»'  and  h2  = 
*^d/ (2+d)  which  make 

— ——  —  (  !  -H  d  -+■  cos  ,t r. Vf  C(  I  d  —  cos  to)) 

(A) 

simple.  C  should  be  chosen  so  that  </»,  —  *  — d  when  a=1(  In  our 
calculation  d=0.55  and  01.942). 

Boundary  Conditions:  conditions  for  the  leading  edge  and 
trailing  edge  stagnation  points  and  the  wall  surface  are 

a  •"  0 .  1  ( </>  “  —  <« ,  *  —  n ,) ;  u  —  0 ,  x  —  £  ]  (11) 

P  — •  0  :-u  —  0,  r  —  4  1 


At  infinity,  B=1 ,  u=U(«p)  and  x=x(£,n=®,  t).  Fortunately,  an 
integral  expression  can  be  found  for  x(5,n  =®,t)  to  simplify  the 
computation  program.  From 


-  '■ ~ — rr  (E.  +  E,  sinV)  -  \EJ\n  **  T'  -  TVcusU  -  d)  (B) 

-  -  -  T_n  \  (X)  J  ■  ? 


we  get 


^  ,  tf,  T  ( </>  ~h  li)  /  2  ]  ~  {  f 

,T,  :n  - ’  —  /: ,!  ca<,{4>  — 

’g  [  (d>»  + 


—  CI  S  (  4>0 


(12) 
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where  E,  =E„+Easin  ol.  <p0=<p(?)  is  already  known.  When  t  is  given, 
and  x(»)  can  be  determined. 

Velocity  distribution  and  displacement  thickness:  when  x 
remains  unchanged,  the  following  can  be  calculated  for  each 
increment 

AC  _  A)' - r,Av,  At,  -  i.AV. 

<iy  (C) 
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Starting  from  the  wall,  £„ «  x,  no=*0,  ^  =C0  -x^Ay  >  >H  =n0  +  (D) 

Vy .  VCn-1  -(x'l)c,-14y-  nn”nn-1+^xe )n-  14y ’  u(W 

u(x,  nAy), .  With  u(x,y),  it  is  not  difficult  to  calculate 

the  displacement  thickness  6,.  . 

III.  Numerical  Method 

Under  the  transformation  (10),  the  basic  equation  (31) 
becomes 

a,  —  'J{0)L " i x'(4>)  +  a^lxiug,  —  2 x.XfU.f  —  *  "...  —  x..uAxtt 

4-  ( x.u t  •+•  xeu.)x.,  —  r^,i«  +  (x.u,  —  xsuj  (  1  3  ) 

•  —  xfxlt/al )]  —  p,u, 

Because  there  is  a  circular  flow  area  on  the  flow  field, 
therefore,  it  is  necessary  to  have  corresponding  countermeasures 
in  the  difference  scheme  and  integration  method. 

1.  Difference  Scheme:  An  implicit  difference  scheme 
(Crank-Nicholson)  is  used  for  partial  derivatives  of  time.  A 
center  difference  scheme  is  used  for  x  ,  x„„,  x  . ,  x  ,  x„ ,  u  , 
uQO,  and  u  Q.  However,  a  tail  wind  difference  scheme  is  used  for 

P  p  ttU 

u  and  u.  terms,  i.e.,  forward  difference  if  the  coefficient  is 
a  p 

positive  and  backward  difference  if  otherwise.  Thus,  the 
difference  equation  can  be  stabilized. 

2.  Because  the  equation  is  non-linear,  its  coefficients  are 
related  to  unknown  functions.  They  must  be  solved  by  iteration. 
We  used  the  Gauss-Sider  iteration  method  to  solve  the  difference 
equation.  In  each  iteration,  n  simultaneous  algebraic  equations 
are  solved  for  a  specific  a  value  (i.e., 


a  specific  x  position). 


n  is  the  number  of  mesh  points  in  the  B-direction.  After  u  and  x 
are  determined,  the  equation  for  the  next  a  is  integrated  to 
complete  the  iteration  step  by  step  from  one  end  to  the  other. 
However,  the  direction  of  integration  of  a  is  changed  after  each 
iteration:  if  it  is  from  o=1  to  a=0 ,  then  the  next  time 
integrate  from  a=0  to  o=1.  Because  the  effect  of  the  circular 
flow  comes  from  both  sides,  this  method  can  quickly  spread  the 
effect  of  boundary  conditions  to  every  point.  Our  experience 
shows  that  this  can  significantly  accelerate  the  rate  of 
convergence . 

3.  Because  the  motion  is  started  abruptly,  it  is  a  singular 
point  at  t=0.  The  proper  selection  of  an  initial  value  when  t  is 
very  small  can  ensure  the  accuracy  of  the  entire  computation. 

The  details  are  discussed  in  reference  [63. 

4.  Determination  of  Point  of  Separation  / 5 75 

The  conditions  for  the  first  occurrences  of  separation 

(singular  point)  are 

n  =  au/ay  =  0  and  D  =  x__x  +  x_  =0  (A) 

€5  nn 

The  deviation  is  in  reference  [63.  However,  D=0  is  no  longer 
true  at  the  point  of  separation  at  later  times.  Therefore,  we 
used  the  direct  condition  (5)  to  determine  the  point  of 
separation.  With  known  x^  and  x^  values  on  a  mesh  point,  the 
intersect  of  x^=0  and  x^=0  can  be  determined  by  intrapolation. 

There  are  usually  two  intersects,  corresponding  to  the  positive 
and  negative  values  of  D.  The  intersect  corresponding  to  D<0  is 
the  true  point  of  separation.  As  shown  in  Figure  5,  there  is 
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x,  —  0, 

another  point  where  x^=0  and  on  the  inside  of  the  enclosed 
curve  (which  corresponds  to  points  outside  the  boundary  layer  as 
discussed  in  3).  This  point  is  the  minimum  of  x(£»n»t).  Thus,  D> 


5.  Choice  of  Mesh  Number:  N_x  N  x  At. 

5  n 

Doubling  or  will  nearly  double  the  computation  time, 
primarily  affects  the  resolution  in  the  5  direction,  which  is 
mainly  determined  by  the  extent  of  pressure  variation.  is 

primarily  determined  by  velocity  distribution  and  by  t  .  For 
instance,  the  circular  flow  region  extends  in  the  n  direction  at 
large  attack  angles.  Therefore,  should  be  increased.  The 
final  and  At  are  determined  by  trial  runs  (i.e.,  doubling 

the  meshes  without  affecting  the  result  significantly). 


l  i 


IV .  Results  and  Discussion 

1.  Trailing  Edge  Separation  -  Leading  Edge  Separation. 

Computations  are  made  at  three  angles  of  attack  (with  the 
semi-axis  ratio  of  the  ellipse  b/a  =  0.5). 

At  zero  angle  of  attack,  1=0,  separation  occurs  at  the 
trailing  edge  at  x  =  2x '  /(a+b)  =  2.85  (the  semi-circumference  is 

S 

3.23).  Several  transformations  a(5,t)  and  B(€,n»t)  were 
attempted  and  results  showed  that  the  one  recommended  in  the 
second  section  is  better.  It  does  not  add  too  much  computing 
time  but  can  maintain  better  accuracy.  The  movement  of  the  point 
of  separation  with  time,  xg~ts,  as  well  as  a  comparison  of  and 
at  separation  and  their  corresponding  known  solutions  are 


shown  in  Figure  1.  Results  obtained  using  two  different  meshes, 
65x17x0.05  and  73x33x0.0125,  are  in  good  agreement.  The  time  at 
which  separation  occurs  is  tg-1.8. 

At  a  medium  angle  of  attack  of  a=22.92°,  separation  still 
occurs  in  the  trailing  edge,  as  shown  in  Figure  2.  Another 
transformation  was  attempted  when  the  angle  of  attack  exists  to 
make  the  mesh  points  closer  near  the  apex  of  the  semi-axis. 


where  C  and  <*„  are  determined  by  a=0  at  <p=-a  and  a-1  at 
However,  the  result  is  not  as  good  as  that  obtained  using 
equation  (10).  At  a  medium  angle  of  attack,  the  pressure 
gradient  near  the  leading  edge  is  not  sufficient  to  cause 
separation.  But,  the  trailing  edge  is  not  sufficient  to  cause 
separation.  But,  the  trailing  edge  stagnation  point  is  shifted 
to  the  flat  position  of  the  ellipse  where  separation  shows  up 
very  late  at  tg=*3.8  because  the  pressure  variation  is  slow.  The 
mesh  remains  to  be  65x17x0.05.  After  tightening  the  mesh  in  the  & 
-direction  by  one  fold,  xg,tg  and  ug  do  not  change  by  one 
percent.  As  the  trailing  edge  singularity  is  developing  (i.e., 
grad  x  is  decreasing  but  not  yet  zero),  the  singularity  near  the 
leading  edge  (where  pressure  varies  violently)  is  also 
developing,  but  much  slower. 

The  separation  at  a  large  angle  of  attack  of  a=50° ,  just  as 
expected,  occurs  at  the  leading  edge  where  pressure  fluctuates 
violently.  The  ratio  of  the  distance  from  the  leading  edge 


stagnation  point  to  the  separation  point,  x(8t)+x  ,  to  the  semi- 

circumference  is  0.32.  Consequently,  the  time  of  separation  also 

arrives  much  earlier;  t  -2.1.  In  this  case,  the  backflow  area 

s 

and  separation  point  are  far  away  from  the  wall.  Therefore,  the 
mesh  is  increased  to  65x65x0.05  (not  too  different  from  the 
result  obtained  with  the  65x33x0.05  mesh).  The  calculated 
results  are  shown  in  Figure  3.  As  the  leading  edge  singularity 
(i.e.,  decreasing  grad  x)  is  developing,  there  is  trailing  edge 
singularity  underdevelopment. 

Summarizing  the  three  cases  above,  separation  always  occurs 

in  the  rear  half  when  &  is  not  large.  While  after  a  exceeds  a 

certain  critical  value  ac>  separation  will  occur  in  the  front 

half.  In  this  case,  tg  will  also  decrease.  Therefore,  the 

variation  of  t  with  the  angle  of  attack  is  as  the  t  ~  &  curve 
s  s 

shown  in  Figure  4.  The  jump  of  the  value  of  t  at  o=o  is  caused 

5>  C 

by  the  "jump"  of  the  point  of  separation  from  the  rear  of  the 
cylinder  to  the  front.  Similar  "jump"  phenomena  also  occur  with 

[ii 

the  steady  separation  of  an  ellipse  .  From  this  study  we  know 

that  singularities  can  develop  in  front  and  rear  of  the  body  due 

to  an  inverse  pressure  gradient  (decreasing  grad  x  with 

increasing  y-axis).  Where  separation  occurs  depends  on  where  the  / 5 76 

singularity  develops  faster.  Singularity  does  not  occur  in  one 

place  and  this  singularity  will  "jump"  from  the  rear  to  the 

front. 

2.  How  to  Improve  Accuracy 

The  special  feature  of  the  Lagrangian  coordinate  is  that  the 


mesh  on  the  (5,n)  plane  varies  with  time.  The  mesh  near  the 
leading  edge  stagnation  point  becomes  larger  as  the  masses  move 
downstream.  Therefore,  the  mesh  should  be  closer.  Next,  in 
order  to  improve  resolution,  the  mesh  should  also  be  closer  near 
the  apex  of  the  a  semi-axis  where  the  pressure  varies  violently. 
However,  too  many  mesh  points  will  lead  to  excessive  computing 
time.  If  the  mesh  is  only  tightened  locally,  not  only  the 
difference  scheme  is  complicated  but  also  the  accuracy  is 
affected  due  to  the  fact  that  the  mesh  is  non-uniform  (for 
instance  the  center  difference  is  no  longer  accurate  to  the 
second  order).  We  solved  the  problem  by  introducing  the 
transformation  (10)  and  an  elliptical  coordinate  (8)  to  make  a 
uniform  mesh  on  the  (a, 8)  plane  to  become  the  non-uniform  mesh  on 
the  (C,n)  plane  -  tightened  at  the  semi-axis  apex  and  the  leading 
edge  stagnation  point.  In  this  case,  although  the  equations  are 
complicated,  yet  logarithmic  calculations  will  not  bring  about 
any  difficulty  in  substance. 


Figure  2  5*22.9° 

o  65X17X0.05  A65X33X0.05 


Figure  3a  5=50° 
o  65X65X0.05  X  33X33X0.1 
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Figure  3b  <**50°,  65X65X0.05 
o  t*2. 1  X  t-2.7 

3.  Variable  Boundary  Problem  -  Steady  Boundary  Problem  / 5 7 7 

In  the  Euler  method,  this  is  a  variable  boundary  problem  on 
the  (x,y)  plane  plane  which  should  be  integrated  from  x=0  to  x=x 

s 

and  x  decreases  with  time.  This  problem  is  very  complex.  In 
s 

the  Lagrangian  method,  this  is  also  a  variable  boundary  problem 
on  the  (£,n)  plane.  However,  it  can  be  treated  as  a  steady 
boundary  problem,  still  integrating  from  a=0  to  a=1.  The  contour 
of  the  solution  x(£,n>t)  after  separation  is  shown  in  Figure  5. 

The  solution  x(£,n,t)  itself  is  normal  on  the  (5,n)  plane.  It 
corresponds  to  finding  y  by  integrating  (4)  along  the  curve 
x*constant  on  the  physical  plane  (x,y).  However,  it  is  not 
possible  to  enter  the  enclosed  curve  x=xg  starting  from  the  solid 
wall,  i.e.,  points  inside  the  curve  are  outside  the  boundary 


layer  (which  is  infinity  on  the  scale  of  V  Re) .  The  x  values  on 
these  points,  in  some  sense,  can  be  considered  as  the  analytical 
extension  of  the  x  values  of  external  points. 

4.  Intensity  of  Separated  Vortex 

At  a  specific  time  t  ,  points  inside  the  singular  curve  x(£, 

,t0)=x  correspond  to  separated  points.  However,  at  a  separated 

singular  point  (50,n0)  9u/ay=0,  u  ^0  and  n  +0  (not  on  the  wall). 

s  s  s  s 

In  this  case,  below  it  we  have  (x=x  ,  n<ne)  8u/ay<0. 

Furthermore,  the  y  value  calculated  according  to  equation  (4)  is 

finite.  Therefore,  this  part  is  not  separated.  Furthermore, 

only  the  part  above  it  is  separated  (x=x  ,  n>n  and  au/ay>0). 

s  s 

According  to  the  M-R-S  condition,  in  the  coordinate  in  motion  at 
(w'^u'-u'),  we  have  w'=0  and  aw'/ay=0.  In  addition,  the  vortex 
flux  separated  from  the  boundary  layer  to  enter  the  external  flow 
is  ^ygw '  aw  * /ay  dy  =  w  ,a  =%(u  '-u '  )a  .  The  intensity  of  the 
separated  vortex  can  thus  be  determined.  It  provides  an 
important  parameter  for  the  next  step  to  take  the  inverse  effect 
of  separated  vortex  on  the  pressure  distribution  into 
consideration.  After  each  step  At,  it  is  possible  to  calculate 
the  newly  increased  separated  vortex  intensity.  (How  to  simulate 
separated  vortex  and  its  movement  is  another  topic.)  Based  on 
the  separated  vortex,  we  can  calculate  a  new  pressure 
distribution.  The  new  pressure  distribution  can  be  used  in  the 
numerical  integration  in  the  next  step  At. 

This  work  was  performed  under  the  guidance  of  Professor  Shen 
Shenfu.  The  author  wishes  to  thank  Dr.  Van  Dommelen  for  his 
unselfish  assistance  during  the  course  of  this  work. 
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Abstract 

.'..i  Fran?’,  in  c  ordinates  we~e  introduced  to  investigate  the  unsteady  boundary  haver 
r"*v.-r»<!o?  of  elliptic  cylinder  a*  angle  of  attack.  The  •  forward-tnu\  mi-nt  of  the  se- 
"r.-a.ti.on  point  with  time  was  e.V«  uhited.  As  well  as  the  first  oeeur.-n .f  separation. 
’*  — "■  rs’so  that  *he  ceparj.*  -vjj  pattern  would  change,  as  exper'-d,  from  ->ar- 

vcpirr.fina  to  fore-separation  as  *he  angle  of  attack  increased. 

In  reject  to  the  mrr.ervr.’  methodogy,  coordinate  transform..;  has  w-re  iritro- 
'’w>'  V>  7V°?  t he  noTu:n;,,''rm  mesh  points  in  the  physical  plane  to  the  unit,  -n  mesh 
~  ■'  *■  ’r.  the  a  txiallir.-y  p  ■»’■<*  s’mplify  the  finite  difference  f  .ra  u’n.  ami  ‘.lie 
moving-boundary  problem  was  treated  in  a  fixed- Viurda*;>  -r. 
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